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Recently observed anomalies in the Bg J/4"t> decay and the like-sign dimuon asymmetry 
hint at possible new physics (NP) in the Bs — B^ mixing. We parameterize the NP with four model- 
independent quantities: the magnitudes and phases of the dispersive part M12 and the absorptive 
part of the NP contribution to the effective Hamiltonian. We constrain these parameters using 
the four observables AMs, AFs, the mixing phase jii^^"^ , and A^j. Our quantitative fit indicates 
that the NP should contribute a significant dispersive as well as absorptive part. In fact, models 
that do not contribute a new absorptive part are disfavored at more than 99% confidence level. 
We extend this formalism to include CPT violation, and show that CPT violation by itself, or 
even in presence of CPT-conserving new physics without an absorptive part, helps only marginally 
in the simultaneous resolution of these anomalies. The NP absorptive contribution to Bs — Bs 
mixing therefore seems to be essential, and would imply a large branching fraction for channels like 

Bs^T+T-. 

PACS numbers; 11.30.Er, 14.40.Nd 
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I. INTRODUCTION 

The Cabibbo-Kobayashi-Maskawa (CKM) paradigm of 
quark mixing in the standard model (SM) is yet to be 
accurately tested in the Bs — Bs sector, and it is quite 
possible that the NP can affect the Bs — Bg system while 
keeping the Bd — Bd system untouched. Indeed, for most 
of the flavor-dependent NP models, the couplings rele- 
vant for the second and third generations of SM fermions 
are much less constrained than those for the first gener- 
ation fermions, alloviring the NP to play a significant role 
in the Bg — Bs mixing, in principle. 

Over the last few years, the Tevatron experiments CDF 
and D0 , and to a smaller extent the B factories Belle 
and BaBar, have provided a lot of data on the Bs me- 
son, most of which are consistent with the SM. There are 
some measurements, though, which show a significant 
deviation from the SM expectations, and hence point to- 
wards new physics (NP). The major ones among these 
are the following, (i) Measurements in the decay mode 

Bs — > J/i^ip yield a large CP-violating phase /3s ^'^'^ [H- 
In addition, though the difference AF^ between the decay 
widths of the mass eigenstates measured in this decay is 
consistent with the SM, it allows AFg values that are al- 
most twice the SM prediction, and also opposite in sign 
. (ii) The like-sign dimuon asymmetry in the com- 
bined B data at D0 0] is almost 4a away from the SM 
expectation. 

The resolutions of the above anomalies, separately or 
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simultaneously, have been discussed in the context of spe- 
cific NP models: a scalar leptoquark model 0,13' mod- 
els with an extra flavor-changing neutral gauge boson 
Z' or R-parity violating supersymmetry Q , two-Higgs 
doublet model [1, Q , models with a fourth generation of 
fermions [13, supersymmetric grand unified models 
p^ . supersymmetric models with split sfermion genera- 
tions [1^ or models with a very light spin-1 particle 
Possible four-fermion effective interactions that are con- 
sistent with the data have been analyzed by [lEj and the 
results are consistent with Q . Similar studies, based on 
the minimal flavor violating (MFV) models , and the 
Randall-Sundrum model [l7[, have been carried out. 

In this paper, we try to determine, in a model- 
independent way, which kind of NP would be able to 
account for both the above anomalies simultaneously. 
We take a somewhat different approach than the refer- 
ences cited above. Rather than confining ourselves to 
specific models, we assume that the NP responsible for 
the anomalies contributes entirely through the Bs — Bs 
mixing, and parameterize it in a model independent man- 
ner through the effective Hamiltonian for the Bs — Bs 
mixing. This effective Hamiltonian "H is a 2 x 2 matrix 
in the fiavor basis, and the relevant NP contribution ap- 
pears in its off-diagonal elements. The NP can then be 
parameterized by using four parameters: the magnitudes 
and phases of the dispersive part and the absorptive part 
of the NP contribution to H. A "scatter-plot" analysis 
that constrained these four new parameters using only 
A'li has been carried out in [l7|. We perform a fit 
to the Bs — Bs mixing observables and obtain a quanti- 
tative measure for which kind of NP is preferred by the 
data. This would lead us to shortlist specific NP models 
that have the desired properties, which can give testable 
predictions for other experiments. It is found that the 
NP needs to contribute to both the dispersive as well as 



2 



absorptive part of the Hamiltonian in order to avoid any 
tension with the data. 

We also extend our framework to include possible CPT 
violation in the Bs — Bg mixing, parameterized through 
the difference in diagonal elements of % . The motivation 
is to check if this can obviate the need for an absorptive 
contribution from the NP. Such an analysis to constrain 
CPT and Lorcntz violating parameters was carried out 
in [l^. However they have used only and not pi^^"^ 
in their analysis, and their parameters are only indirectly 
connected to the elements of T-L. We try to account for 
the two anomalies above with only CPT violation as the 
source of NP, and with a combination of CPT violation 
and the NP contribution to the off-diagonal elements of 
T-L. As we will show, nothing improves the fit significantly 
from the SM unless there is a nonzero absorptive part in 
the Bs — Bs mixing amplitude. 

The paper is organized as follows. In Sec. we in- 
troduce our formalism for the four NP parameters. In 
Sec. IIII[ we summarize the experimental measurements 
and theoretical predictions for the observables relevant 
for Bs — Bs mixing. In Sec. IIV| we present the results 
of our fits, and their implications for NP models are dis- 
cussed in Sec. |Vl In Sec. IVIl we introduce the formalism 
for introducing CPT violation and in Sec. lVIll we explore 
the extent to which it can help resolving the anomalies. 
Sec. I Villi summarizes our results and concludes. 

II. THE EFFECTIVE HAMILTONIAN 



Arg(r?f )■ We take the phases Arg(MNP) and Arg(rJ'2P) 
to lie in the range Q-2tt. 

In a large class of models, including the Minimal Flavor 
Violation (MFV) models, the NP contribution has no 
absorptive part, i.e. T12 = P™. This is true for a lot 
of non-MFV models too. This occurs when NP does not 
give rise to any new intermediate light states to which 
Bs or Bs can decay. For such models, Eq. ^ implies 
that AF^ < Ars(SM) « 2|rf^|, i.e. the value of AF^ is 
always less than its SM prediction [l^. In such models, 
the NP is parameterized by only two parameters: |M;f^P| 
and Arg(Af}^P). An analysis restricted to this class of 
models was performed in [20|. 

However there exists a complementary class of viable 
models where the NP contributes to F12 substantially. 
These include models with Icptoquarks, R-parity violat- 
ing supersymmetry, a light gauge boson, etc. It has been 
pointed out in Q that such a nonzero absorptive part 
that arises naturally in these class of models can en- 
hance AFg significantly above its SM value, contrary to 
the popular expectations based on [l^. One notes that 
a new absorptive part in the mixing amplitude necessar- 
ily means new final states that can be accessed by both 
Bs and Bs- The data from the direct measurements of 
branching ratios is extremely restrictive [Tsj , apart from 
that for a few final states like Bs — >■ r+r" Q. As we 
shall see later in this paper, such models are favored by 
the Bs — Bs mixing data. The importance of t+t^ final 
states from Bd and Bs decays has also been pointed out 
in M- 



The evolution of a Bg — Bs state can be described by 
the effective Hamiltonian 
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in the flavor basis, where Mij and F^ are its dispersive 
and absorptive parts, respectively. When CPT is con- 
served, Mil = M22 and Fn = F22. The eigenstates of 
this Hamiltonian are Bsh and Bsl , with masses Msh and 
MsL respectively, and decay widths Tsh and Fsl respec- 
tively. The difference in the masses and decay widths can 
be written in terms of the elements of the Hamiltonian 
as 

AAf, = MsH - MsL ~ 2\Mi2\ , 
AF, = TsL - TsH ~ 2IF12I cos[Arg(-Afi2/Fi2)] .(2) 

The above expressions are valid as long as AFg ^ Ms, 
which is indeed the case here. 

Since CPT is conserved, the effect of NP can be felt 
only through the off-diagonal elements of "H. We separate 
the SM and NP contributions to these terms via 



M12 = Mf^ + M^i" 

_ -pSM I T^NP 
-1-12 — -1-12 + -1 12 • 



(3) 



The NP can then be completely parameterized in terms 
of four real numbers: |Mf/|, Arg(Mf/), \r^^\ and 



III. THE MEASUREMENTS 

The Bs — Bs oscillation and CP violation therein can 
be quantified by four observables, viz. the mass difference 
AMs, the decay width difference AF^, the CP-violating 
phase pi^^'^ , and the semileptonic asymmetry afj. 

The mass difference is measured to be 



AMs = (17.77 ± O.IO ± 0.07) ps'^ , 
which is consistent with the SM expectation [22 
AAfg(SM) (17.3 ± 2.6) ps^^ . 



(4) 



(5) 



However measurements in the Bs — > J/iJj4' decay mode 
show a hint of some deviation from the SM. The CP- 
violating phase jS's^'^'^ in this decay is 
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whose average value measured at the Tevatron experi- 
ments m is 



(0.4itO;}Du(i.i6lS:}^3) 



In the SM, 



^J/^0(SM) = Arg ( - 
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0.019 ±0.001 



(7) 
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Thus, the measured value of jSg^^'^ is more than 2(t away 
from the SM expectation. On the other hand, the differ- 
ence in the decay widths of the mass eigenstates and 
is measured to be [l| 



which is about 1.7a away from the SM prediction 



Ar, - ±(o.i54t0:°57^) ps-i , 



while the SM expectation is [22 1 

Ar^(SM) = (0.087 ±0.021) ps^^ 



(9) 



(10) 



The measurement is consistent with the SM expectation 
to ~ la, however it allows for AF^ values that are al- 
most twice the SM prediction. Note that the sign of AF^ 
is undetermined experimentally and this gives us more 
room to play with the NP parameters. 

CDF has recently announced its new results, based on 
5.2 fb-i of data ^ 

\Ar,\ = (0.075 ± 0.035 ± 0.010) ps^S 

f^J/i^<t> ^ (0.02 - 0.52) U (1.08 - 1.55) (11) 

to 68% C.L.. While we note that the results are consis- 
tent with the SM, the final Tevatron averages are still 
awaited. Therefore, we use the values in Eq. ^ in our 
analysis. 

The other anomalous measurement is the like-sign 
dimuon asymmetry. Averaging the 9.0 fb^^ data of 
D0 g and 1.6 fb-i data of CDF [2|, and adding the 
errors in quadrature and treating them as Gaussian, we 
get 



= -(7.41 ±1.93) X 10" 



(12) 



which differs by more than 3ct from its SM prediction 

A^,(SM) = (-0.23i;5:°^)xlO-3. (13) 

Note that for A^j, CDF has a poorer statistics than 
D0 and therefore the average value is dominated by the 
D0 data. 

Even in the presence of new physics, the SM relation- 
ship holds: 

All = (0-506 ± 0.043)a^i + (0.494 T 0.043X1 , (14) 

where a|i and a^i are the semileptonic asymmetries for the 
Bs — Bs and the Bd — Bd systems, respectively. The for- 
mer is related to the Bs — Bs mixing obscrvablcs through 
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where (fis = Arg(— Afi2/Fi2). The latter is defined anal- 
ogously. The coefficients in Eq. are experimen- 
tally measured, and contain information about AM^jj^-j, 
AFrf(s), and production fractions of Bd and Bs mesons. 
Using afi = -(4.7 ± 4.6) x lO'^ 0, this leads to 



<i(SM) = (2.06 ±0.57) X 10 
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(17) 



The value of afj depends on AAf^i, AF^ and (pd, the pa- 



rameters in the Bd sector anologous to those in Eq. (|15p . 
These parameters depend on the NP in the Bd sector, 
which is independent of the NP parameters in the Bs 
sector that we are considering. We therefore do not con- 
sider the measured values of as a direct constraint, 
but express it in terms of A Aid, ATd, and 0^, whose 
experimental values are taken as inputs. 

In the SM, we have 0s(SM) = 0.0041 ± 0.0007 
Note that if the dominating contribution to Fi2s were 
from a pair of intermediate c quarks, (j)s(SM) would have 
been equal to —2/?/^''^'^. Since the intermediate u — c and 
u — u quark states give comparable contributions to Fi2s, 
we have 0,(SM) ^ -2/3s''^'^'^(SM) [H. 



IV. THE STATISTICAL ANALYSIS 



AM,,AF3,/3//'^'^ and a: 



all = -0.010 ±0.006 , 



(16) 



We perform a fit to the observed quantities 
gi, using the NP parameters 
Arg(MNP), \r^P\ and Arg(rJ^2P). We assume 
all the measurements to be independent for simplicity, 
though the measurements of AF^ and pi^^'^ are some- 
what correlated. The values of all the obscrvablcs and 
their SM values are as given in Sec. IIIII In order to ex- 
press them in terms of M12, -^^2'^, F12 and Ff|^, one has 
to use Eq. ^ in addition. In order to take into account 
the errors on the SM parameters, we add the theoreti- 
cal and experimental errors on our observed quantities 
in quadrature. 

Note that since we have four observable quantities and 
four parameters, it is not surprising that we obtain the 
global minimum value of as x^i„ = when all the 
NP parameters are allowed to vary. The questions we 
address here are (i) what the preferred values of the NP 
parameters are, and (ii) to what confidence level (C.L.) 
a given set of NP parameters (or SM, which is a special 
case of NP with Mf^ = T^:^ = 0) is allowed. The lat- 
ter is obtained assuming all errors to be Gaussian. Here 
we give our results in terms of the goodness-of-fit con- 
tours for the joint estimations of two parameters at a 
time. The [la, la, "ia, 4(t) contours, that are equivalent to 
p-values of (0.3173,0.0455,0.0027,0.0001), or confidence 
levels of (68.27%, 95.45%, 99.73%, 99.99%), correspond to 
^ = (2.295,6.18, 11.83, 19.35), respectively. 

In Fig. [I] we show the la, 2a, 3a, 4a contours in the 
IM12I — Arg(Mi2) plane, where the other NP parame- 
ters are marginalized over. Clearly, we see a preference 
towards nonzero jA/j'^^l as well as nonzero ATg{Mf^) 
values. There are two best-fit points with = 0, 
one at Mf.^ « 6.3exp(2.0 i) ps~^ and the other at 
« 16.2 exp(2.8 i) ps~^, shown with crosses in Fig.[T] 
Actually, each of these crosses is a superimposed double, 
with two values of T^2 i as shown in Fig. ^ The points 
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FIG. 1: The la (red/solid), 2a (green/dashed), 3a 
(blue/dotted) and 4(t (pink/dot-dashed) goodness-of-fit con- 
tours in the \Mi^\ - Arg{Mi^) plane, where the other NP 
parameters are marginahzed over. The best-fit points, with 
= 0, are denoted by crosses. 



correspond to the constructive and destructive interfer- 
ence between the SM and NP amphtudes in order to 
give the measured central values of AMg. The region 
with Mf^ — 0, i.e. the x-axis, is outside the 2cr re- 
gion, indicating that it will be rather difficult to fit the 
current data without some NP contribution to the dis- 
persive part of the Bg — Bg mixing. The contours also 
imply that \M^^\ < 21.1 ps'^ to 3a. 

In Fig. [21 we show the goodness-of-fit contours in the 
|ri2| — Arg(ri2) plane, marginalizing over other two NP 
parameters. As the measurements do not determine the 
sign of APs, for any particular value of lAPsj, we per- 
form the fit for both positive and negative values. 



and keep the minimum of the two. This doubles the 
number of best-fit solutions, and the two best-fit points 
of Fig. [T] now split into four. For |A/^^^| ~ 6.3, the 
solutions are Pjif = 0.18exp(6.0 i) or 0.18exp(5.2 i), 
and for |Af{^^| = 16.2, the corresponding solutions are 
Tf.[ = 0.18exp(0.2 i) or 0.18exp(l.l i) (both Aff/ and 
T^2 are in ps^^, here, and also later where not men- 
tioned explicitly) . Note that there is a reflection symme- 
try about Arg(P5^2^) = tt. Again, a preference for nonzero 
values of |r5^2^| is indicated, though Arg(P5^2^) may van- 
ish. The region with T^^ = 0, i.e. the X-axis, is outside 
the 4(7 allowed region, indicating that NP contribution to 
the absorptive part of the effective Hamiltonian is highly 
favored. The contours also imply that \Ti2\ ^ 0.35 at 
3cr. 
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FIG. 3: The la (red/sohd), 2a (green/dashed), 3cr 
(blue/dotted) and 4(t (pink/dot-dashed) goodness-of-fit con- 
tours in the |Af^''| — \T^2\ plane, where the other NP param- 
eters are marginalized over. The best-fit points, with ~ 0, 
are denoted by crosses. 
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FIG. 2: The la (red/solid), 2a (green/dashed), 3a 
(blue/dotted) and 4cr (pink/dot-dashed) goodness-of-fit con- 
tours in the \rii\ - Arg(r^^f ) plane, where the other NP 
parameters are marginalized over. The best-fit points, with 
= 0, are denoted by crosses. 



Fig. [3] displays the contours in the | Af^^ | — |P5^2^ | plane, 
and the two NP phases are marginalized over. Not only 
does it show a preference for nonzero values of M^.^ and 
Ti2, but the A/^^ = axis is outside the 2a allowed 
region and the = ^^^^ outside the 4a allowed 
region. The best fit points are again superimposed dou- 
bles, whose values can be read off from the discussion 
above. The origin in this figure is the SM, which has 
XsM ~ 25.85, and lies even outside the 4cr allowed re- 
gion. This dramatically quantifies the failure of the SM 
to accommodate the current data. The reason is evident 
from eqs. (0) and ([15]); while B^ -> J/tPcj) prefers /j/^'^'^ 
close to 7r/8 or 37r/8, with a probability minimum near 

pj/'^'^ Ri 7r/4, the measurement of A^^, and hence that of 

a|i, prefers large tan^^, forcing pi^'^'^ close to 7r/4. This 
creates the tension between these two measurements. 

Fig. [3] also tells us that the models for which T^2 = 0; 
like R-parity conserving supersymmetry, universal extra 
dimension, and extra scalars, fermions, or gauge bosons, 
cannot bring the tension down even to the 4(7 range, un- 
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less the data moves towards the SM expectations (and 
unless the new bosons are flavor-changing so as to gener- 
ate a nonzero r5^2^). The best fit point with Ffj^ = has 
= 20.75 and corresponds to Mf^ = 3.72exp(1.68 i). 
This is further emphasized in fig. |4l which shows the 5ct 
[p value of 10~^, = 27) contour for those NP mod- 
els where r5^2^ is set to vanish (within the closed contour 
above and under the open contour below). 
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FIG. 4: The Scr goodness-of-fit contour in the |Mf^^| - 
Arg(Afi'^^) plane, when V12 = 0, i.e. NP does not contribute 
to the absorptive part of the effective Hamiltonian. There are 
no points that are allowed to within 4cr. The best-fit point, 
with = 20.75, is denoted by a cross. 

One may question the optimistic SM uncertainty for 
AFs as quoted in eq. pU)) . However, this has an almost 
negligible effect. For example, if we increase the uncer- 
tainty by 50%, neither the best fit points nor the con- 
fidence levels change significantly. The best fit point in 
fig. Ulhas a x^ minimum of 20.58 instead of 20.75. The 
reason is the large deviation of from its SM value, 
to explain which we need a significant enhancement in 
tan <h^. 



V. PREFERRED NP MODELS 

From the results and discussion in the previous section, 
it appears that: 

(i) The SM by itself is strongly disfavored. Either Mf^ 
or should be nonzero. 

(ii) Mf^ ^ but Fjlf = is also not allowed at 4cr, but 
the fit is marginally better than the SM. 

(iii) The hypothetical case where Fjif 7^ but A/f/ = 
is also disfavored to more than 2a. (This is a rather 
natural condition, since any interaction that contributes 
to F5^2^ will necessarily contribute to Af^^.) 

Most of the NP models can contribute significantly to 
M^2 ■ Leading examples are the MFV models like min- 
imal supersymmctry, universal extra dimensions, little- 
Higgs with T-parity, etc. Non-MFV models like a fourth 
chiral generation, supersymmctry with R-parity viola- 



tion, two-Higgs doublet models, models with extra Z', 
etc. can also contribute significantly to ■ 

The NP models that can contribute significantly to 
T^2 ^ however, arc rather rare. This is because the NP 
contribution to the absorptive part needs light particles 
in the final state, and there are strong limits on the de- 
cays of Bg to most of the possible light final state par- 
ticles. One of the few exceptions is the mode r+r", on 
which there is no available bound at this moment. Thus, 
the NP that contributes to T^2 to do so via the inter- 
action h — > sT+r^, but without affecting related decays 
like h — se'^e~ or 6 — > s^^/x" [Tsj . This can be achieved 
only in a limited subset of models, for example those with 
second and third generation scalar leptoquarks, R-parity 
violating supersymmctry Q, or extra Z' bosons It 
turns out that the former can provide enough contribu- 
tion to F5^2^ to increase AFg up to its current experimen- 
tal upper bound 0, @ . The amount of NP required for 
this is consistent with the difference between the decay 
widths of Bd and B, mesons (r,,/F<j - 1 = (3.6 ± 1.8)% 
) , and the recent measurement of the branching ratio 
of B+ K+T+T-, which is less than 3.3 x lO'^ at 90% 
C.L. H^. 

One should note here that if the D0 results on the 
dimuon charge asymmetry survive the test of time, it 
will be a clear indication of the presence of a nonzero 
r5^2^. Such models are also favoured from the CDF and 
D0 combined result on the allowed contours for pi^^'^ 
and Fg, but we need to wait for the final Tevatron aver- 
age. 



VI. CPT VIOLATION: THE FORMALISM 

The analysis till now is valid only if we assume CPT- 
invariance. However, the CPT symmetry may be violated 
in theories that break Lorentz invariance [27[ . Indeed for 
local field theories, CPT violation requires Lorentz vi- 
olation [23. (This need not be true for nonlocal field 
theories as well as for theories with noncommutative 
space-time geometry, see [1^.) In general, CPT violation 
should result in differences in masses and decay widths 
between particle-antiparticles pairs. However it may be 
easier to identify even through oscillation experiments, 
which typically arc sensitive to an interference between 
the CPT-conserving and CPT-violating interactions. 

While CPT violation in the K system is severely con- 
strained through the mass difference between the neutral 
kaons [s^l, the bounds on the CPT violating parame- 
ters in the Bd and Bs systems arc rather weak. In fact, 
the bounds for the Bd sector are about three orders of 
magnitude weaker than those for the K sector [3l[. The 
bounds on Lorentz-violating parameters using the data 
on B mesons can be found in |18j and references therein. 
Here we use a model-independent parameterization, like 
the one earlier followed in [32j and recently used by two 
of us (ssj , and determine the preferred parameter space 
using the data on Bs — Bg oscillations. Unlike [l^, we 
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take both A^^ and pi^^"^ data into account. 

One should note that as a new physics option, CPT vi- 
olation is not exactly at the same footing as the models 
mentioned before. However in the language of the effec- 
tive Hamiltonian the CPT violation manifests itself 
naturally through in a difference between the diagonal 
elements oi H. It is therefore interesting to see if the 
constraints on the NP coming from AM^ and AF^ can 
be relaxed at all with these additional degrees of freedom. 
A posteriori, we will justify the discussion on CPT viola- 
tion by showing that if the new physics indeed turns out 
to be without an absorptive part, CPT violation might 
help to explain the Bs — Bg mixing data, albeit only 
marginally. 

The CPT violation manifests itself in the effective 
Hamiltonian through the difference in the diagonal el- 
ements. We write the effective Hamiltonian in eq. ([T]) 



n = 



Mo - iPo - 6' Mi2 2 
M^2 - 1^2 Mo - fPo + S' 



(18) 



and define the dimensionless CPT- violating complex pa- 
rameter S as 



25' 



V H12H21 \J H12H21 



(19) 



where Hij = Mij 



-P 



The eigenvalues of H are 



Mo - ^Po ) ± ayHi2 , 



(20) 



where a = \/ H21/H12 and y = \J\ + (5^/4 . The corre- 
sponding mass eigenstatcs are 



with |pi 



\Bsl) - V2 \Bs) - q2 \Bs) , 

k2| 



(21) 




(22) 



Clearly, CPT invariance corresponds to ?/i = 1]2 ■ 

Let us now determine the dependence of our four ob- 
servables on the CPT-violating parameters. The differ- 
ences in masses and decay widths of the eigenstates are 
related to the difference in eigenvalues as 



Ai - A2 = AM + -AT , 



(23) 



where Ai and A2 arc ordered such that Rc(Ai — A2) > 0. 
From Eq. (|20|) . 

AM Ml - M2 = 2Rc{ayHi2) , (24) 
AF = F2 - Fi = 41m{ayHi2) • (25) 



Since IF12I ^ IM12I, we can write 



aHi2 = I A/12 1 
« I A/12 1 



1 - 



1 |ri2p 



2 VAfi2 



iRc 



ri2 
A/12 



Then Eqs. ^ and ^ yield 



AM 
AF 



lA/i 



IM- 



121 



2Rc(?/) + Im(y)Re 
4Im(y) - 2Re(y)Re 



Fi 



M12 
"F12 



M- 



12 



(26) 

(27) 
• (28) 



The dependence on the CPT-violating parameter S ap- 
pears entirely through y. 

Let us pause here for a moment and find what the 
above two equations tell us about the allowed parame- 
ter space. Let us first focus on the best constraint, AMg, 
and work in the limit where F12/A/12 is negligible. IA/12I, 
and hence Mf^, can be arbitrarily large, as Re(y) can 
be made arbitrarily small by an appropriate choice of 
6. Similarly, Rc(j/) can be quite large (albeit compatible 
with other constraints) as long as there is a near-perfect 
cancellation between the SM and NP mixing amplitudes, 
making \Mi2\ small. However, the smallness of AT/ AM 
constrains Im(j/) /Re(y) to be small, thus indicating that 
y is almost real. Since y = \/l + (5^/4, this implies that 
6^ is almost real and Re((5^) > —4. Therefore, one would 
expect that S is either almost real, or it is almost imagi- 
nary, but with |Im((5)| < 2. 

Now let us consider the CP-violating observables (Sj^^'^ 
and a||. The effective value of the former may be ob- 
tained in the presence of CPT violation by considering 
the decay rates of Bg and Bg to a final CP eigenstate 
fcp as [H: 



T{Bs{t) ^ fcp) 
r(B4t) ^ fcp) 



\Af\'[\U{tr + \^f,\'\f-it)\' + 
2Re{^fJ.{t)fl{t))] , (29) 

i^/i2r 



'■[\f-{t)\' + \O.W+it)\' + 

2Re(e/J+ (0/1(0)] , (30) 



with 



^ . (31) 

V2 



Here Af and Af are the amplitudes for the processes 
Bs ^ fcp and Bg — >■ fcp, respectively. The time evolu- 
tions are given by 



1 



1+OJ 

1 



/+(*) = -— (e-'^^*-Hwe-'^^*), 
1 + U) 



Uit) 



1 



1+OJ 



(32) 
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The final state in Bg — >■ J/ ■00 is not a CP cigenstate, but 
a combination of CP-even and CP-odd final states, which 
may be separated using angular distributions. With the 
transversity angle distribution [s^], the time-dependent 
decay rate to the CP-even state is given by the coefficient 
of (1+cos^ 9), while the time-dependent decay rate to the 
CP-odd state is given by the coefficient of sin^ 6. 

The value of effective pi^^'^ in this process is deter- 
mined by writing the time evolutions and ([5(1)) in 
the form 

T(B,{t) ^ fcp) = ci cosh(Ar,t/2) -f- C2 sinh(Ar,i/2) + 
C3 cos{AMst) + C4 sin{AMst) , (33) 

T(Bs{t) fcp) = ci cosh(Ar,t/2) + C2 sinh(Ar3t/2) -t- 
C3 cos(AA/^i) + C4 sin(AM^t) . (34) 



The direct CP violation in Eg is negligible; i.e. 

1. Also, \Ti2/Mi2\ < 1, so that in the absence 
of CPT violation, |?7i| = \ri2\ ~ 1. Then in terms of 
= = = aAf/Af, one can write 

where i]cp is the CP eigenvalue of fcp- 

When CPT is violated, the effective phases pi^^'^ and 
pi^^'^ measured through Bs{i) and Bs{t) decays, respec- 
tively, will turn out to be different. Indeed, the difference 
between these effective phases will be a clean signal of 
CPT violation. 



sin(2/3//'^^) 
sin(2^//'^^) 



-riCP-, 



2[-Im(w) - Re(^/JIm(c.) + Im(OJ + Im(e/jRe(w)] 



-■ncp 



2h|0,plm(c.) 



• 2Re(^/J - 2Re(0JRe(w) 
+ Re(e/JIm(^)+Im(^/,)- 



- 2Im(^;JIm(w)] 
Im(e/jRe(^)] 



[2 + |^/.P(1 + I^P) - 2Rc(0J + 2Re(0,)Re(a.) - 2Im(C;JIm(c.)] 



(36) 
(37) 



Though the analysis of the Bs and Bg modes needs to be 
performed separately, here we assume identical detection 
and tagging efficiencies for both, and use the average of 
Eq. ^ and Eq. ^ for our fit. 

The semileptonic CP asymmetry a|j is measured 
through the "wrong-sign" lepton signal: 

^ T{Bs{t)^^i+x)~T{BS)^^^-x) _ ^gg^ 



T{Bs{t) -> ^i+X) + T{Bg{t) ^ X) 



Here, 



T(Bs{t)^ ^i-X) = \inif^A{Bs^ ^i+X)W (39) 
T{Bg{t) ^ ^i+X) = \{f^/r^2)A(Bg ^ ^i+X)\^ ,(40) 



and since \A{Bs M+X)| = \A{Bs M+X)|, 



1 

ii)2r 



l^il 



si ~ 1 I I 19 



1 



(41) 



which is independent of the CPT-violating parameter 
5. That the semileptonic asymmetry does not contain 
a CPT violating term in the leading order was also noted 
earlier 1351. 



VII. CPT VIOLATION: THE STATISTICAL 
ANALYSIS 

In this Section, wc perform a to the observables 

and a|j . Let us first 



AM<t, AFs, the effective phase 



assume that there is no CPT-conserving NP contribution 
coming from M^2 ^^^"^ ^^2^, so that the only relevant NP 
contribution is CPT violating, and is parameterized by 
Rc((5) and Im((5). The allowed parameter space is shown 
in Fig. [5j It turns out that in this case, the value of Xmin 
is « 16.4 (at 5 = 0.008+0.958 i and 5 = -0.024-h0.958 i), 
marginally better than the one obtained in the {T^2 — 
0, 7^ 0) case discussed above in Fig. U) There are 

some, albeit small, regions in the parameter space that 
are allowed to 4(7. However a fit good to 3(t or better is 
still not possible. 

We therefore need to add the CPT-conserving NP to 
the CPT-violating contribution. However we have al- 
ready seen in the preceding section that and 
together are capable of explaining the data by them- 
selves. Therefore the fit using (5, Mf^ as well as T^2 is 
redundant. With six independent parameters and only 
four observables, not only is x^in = guaranteed, but 
no effective limits on CPT-conserving and CPT-violating 
parameters arc generated. 

Wc. therefore, go directly to the possibility where there 
is CPT-conserving NP, but without an absorptive part: 
= 0. We have already observed (Fig. H]) that the 
entire region in the — Arg(M|^^) is outside the 

4(7 region in such a scenario. We would now ask what 
happens if we enhance the two-parameter NP with two 
more CPT violating parameters, viz., Re((5) and Im((5). 
This scenario is interesting because, as we have seen be- 
fore, only very specific kind of NP can contribute to T^2 ^ 
which would be tested severely in near future. In case no 
evidence for the relevant NP is found (e.g. the branching 
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FIG. 5: The 4a goodness-of-fit contours in the Re(5) — lm{S) 
plane, when the only relevant NP contribution is CPT violat- 
ing, parameterized entirely by 5. There are no points that are 
allowed to within 3a. The crosses show the best fit points. 



with X 



16.4. 



ratio of Bg — > r+r" is observed to be the same as its 
SM prediction), the next step would be to check if CPT 
violation, along with the NP contribution through , 
would be able to account for the anomalies. For example, 
one may want to determine and /3s of Eqs. ([36]) and 
(f37| separately and see whether they are different. 

Fig. El shows the situation in the \M^^\ - Arg(Mf^P) 
plane. As compared to Fig. 21 one can see that once we 
marginalize over S, we now have some regions allowed to 
within 4(7 (within the closed contour above and below the 
open contour), but none within 3a. Indeed, Xmin = 14.3 
at = 3.54 exp(5.76 i). This clearly does not improve 
the goodncss-of-fit substantially, indicating that there is 
no good alternative for rf2^. 
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FIG. 6: The 4a goodness-of-fit contours in the |Mi2 | - 
Arg(AfiY) plane, when Ffa^ = 0, i.e. NP does not con- 
tribute to the absorptive part of the effective Hamiltonian. 
The CPT- violating complex parameter 5 has been marginal- 
ized over. There are no points that are allowed to within 3a. 
The cross shows the best fit point, with = 14.3. 



Fig.[2]shows the situation in the complex 5 plane, when 
has been marginalized over. The best-fit point cor- 
responds to S = —0.01 + 1.40 i, which gives Xmin — 14.3 
as mentioned earlier. The CPT conserving point {S = 0) 
lies outside the 4(t region. As expected from the discus- 
sion in Sec. IVIl the allowed values of 5 arc close to the 
Rc(5) or Im((5) axis, with |Im((5)| restricted to 2. One 
observes that the current data allows rather large 1) 
positive values of Im((5) at 4a. 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 
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Re[5] 
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FIG. 7: The 4a goodness-of-fit contours in the Re(5) — lm{S) 
plane, when the complex NP parameter M^2' is marginalized 
over, while r5^2^ has been constrained to vanish. There are 
no points that are allowed to within 3a. The cross shows the 
best fit point, with = 14.3. 



VIII. CONCLUSION 

Any flavor-dependent new physics model can in gen- 
eral affect both mass and width differences in the B^- 
Bs system. It can also affect the CP-violating phase, as 
well as the dimuon asymmetry, which was found by the 
D0 collaboration to have an anomalously large value. 
With these four obscrvablcs, one can constrain the free 
parameters of the new physics model. We have used the 
model independent approach where we consider the ef- 
fective Bs — Bs mixing Hamiltonian T-L and parameterize 
the NP through its contribution to T-L. We quantify the 
goodness-of-fit for the SM and NP parameter values by 
performing a combined x^-fit to all the four measure- 
ments. The tension of the data with the SM is clear by 
the high value of at the SM. Moreover, it is observed 
that we need NP to contribute to the dispersive as well 
as absorptive part of the off-diagonal elements of T-L in 
order for the current data to be explained. The absorp- 
tive contribution, in particular, can be obtained from a 
very limited set of models, which will be severely tested 
in near future. 

We also introduce the possibility of CPT violation by 
adding unequal NP contributions to the diagonal ele- 
ments of "H. We explicitly show how CPT violation might 
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affect the observables, especially dwelling on the effect on 
pJ/^'t' ^ Taken alone, the CPT violation cannot affect the 
dimuon asymmetry, and it can make the fit to the Bs — Bg 
mixing data only marginally better. In combination with 
a CPT conserving NP, it can enhance the allowed param- 
eter space for that NP, however it docs not seem to be 
able to obviate the need of an absorptive contribution 
from NP. 

The data on all the observables considered in this pa- 
per is still relatively preliminary, the deviations from the 
SM are only at about 2-3ct level, and future data may 
either confirm these deviations or expose them as statis- 
tical fluctuations. If the errors and uncertainties shrink 
keeping the central values more or less intact, this will 
mean: 

• The SM is strongly disfavored. Moreover, the rele- 
vant NP should be flavor-dependent, as we do not 
see much deviation in the B^ — Bd sector. 

• The NP models that do not contribute to the ab- 
sorptive amplitude of the Bg — Bg mixing are also 
strongly disfavored if CPT is conserved. The best 
bets are those NP models that provide both dis- 
persive and absorptive amplitudes in the Bg — Bg 
mixing. This also gives rise to new decay channels 
for Bg. For example, one might find the branching 
ratio of Bg — > r+T~ enhanced significantly from its 



SM expectation. 

• Without any CPT-conserving NP, only CPT viola- 
tion is only of marginal help, as it cannot enhance 
the semileptonic asymmetry. Even in combination 
with the CPT-conserving dispersive NP, it cannot 
allow regions in the parameter space to better than 
3cr. 

To summarize, the NP models that contribute an ab- 
sorptive part to Bg — Bg mixing seem to be essential if 
one wants to explain the data on and simulta- 

neously. There is only a limited set of such models, and 
they will be severely tested in near future. In the sce- 
nario that such an absorptive NP contribution is ruled 
out, one may have to resort to CPT violation in order to 
explain the data. A prominent signature of such a CPT 
violation would be a difference in jSg^^'^ and pi^^'^ as 
shown in eqs. (|36|) and (|37| . 
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